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Abstract. In this work, we show the convergence of adaptive lowest-order FEM (AFEM) 
for an elliptic obstacle problem with non-homogeneous Dirichlet data, where the obstacle 
X is restricted only by \ £ H 2 (£l). The adaptive loop is steered by some residual based 
error estimator introduced in Braess, Carstensen & Hoppe (2007) that is extended to 
control oscillations of the Dirichlet data, as well. In the spirit of Cascon et AL. (2008), 
we show that a weighted sum of energy error, estimator, and Dirichlet oscillations satisfies 
a contraction property up to certain vanishing energy contributions. This result extends 
the analysis of Braess, Carstensen & Hoppe (2007) and Page & Praetorius (2010) 
to the case of non-homogeneous Dirichlet data as well as certain non-afHne obstacles and 
introduces some energy estimates to overcome the lack of ncstcdncss of the discrete spaces. 



1. Introduction 



1.1. Comments on prior work. Adaptive finite element methods based on various 
types of a posteriori error estimators are a famous tool in science and engineering and are 
used to deal with a wide range of problems. As far as elliptic boundary value problems are 
concerned, convergence and even quasi-optimality of the adaptive scheme is well understood 
and analyzed, see e.g. [4, 15, 18, 28, 29, 36, 37]. 

In recent years the analysis has been extended and adapted to cover more general applica- 
tions, such as the p-Laplacian [39], mixed methods [12], non-conforming elements [13], and 
obstacle problems. The latter is a classic introductory example to study variational inequal- 
ities which represent a whole class of problems that often arise in physical and economical 
context. One major application is the oscillation of a membrane that must stay above a 
certain obstacle. Other examples are filtration in porous media or the Stefan problem (i.e. 
melting solids). In both of which, non- homogeneous Dirichlet data play an important role. 
Also in the financial world, obstacle problems arise, e.g. in the valuation of the American put 
option [33], where one has to deal with various non-affine obstacles. For a broader under- 
standing of these problems, we refer to [20] and the references therein. The great applicability 
in many scientific areas thus make numerical analysis and mathematical understanding of 
the obstacle problem both, interesting and important. As far as a posteriori error analysis is 
concerned, we refer to [5, 7, 8, 16, 25, 30, 38]. Convergence of an adaptive method for elliptic 
obstacle problems with globally affine obstacle was proven in [9, 32]. Both of these works, 
however, considered homogeneous Dirichlet boundary data and affine obstacles only. 
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1.2. Contributions of current work. We treat the case of a general obstacle x £ H 2 (Q). 
By a simple transformation and allowing non- homogeneous Dirichlet data (Prop. 4), this 
can, however, be reduced to the case of a constant zero-obstacle. Since our analysis works for 
general globally affine obstacles, even without the reduction step, we consider affine obstacles 
and non-homogeneous Dirichlet data in the following. We follow the ideas from [32], i.e. 
adaptive Pl-FEM for some elliptic obstacle problem with globally affine obstacle. Contrary 
to [32] and [9], however, we allow non-homogeneous Dirichlet boundary data g G H l (T), 
which are approximated by gt via nodal interpolation within each step of the adaptive loop. 
In contrast to the aforementioned works, we thus do not have nestedness of the discrete 
ansatz sets, which is a crucial ingredient of the prior convergence proofs. In the spirit of [15] 
and in analogy to [32], we show that our adaptive algorithm, steered by some estimator Q£, 
guarantees that the combined error quantity 

(1) A e := J(U e ) - J(u e ) +lQ* + \ apx 2 
is a contraction up to some vanishing perturbations at — > 0, i.e. 

(2) A m < kA £ + a e , 

with 0<7, At < 1, A > 0, and at > 0. The data oscillations on the Dirichlet boundary are 
controlled by the term apx £ , and the quantity ut denotes the continuous solution subject to 
discrete boundary data which is introduced to circumvent the lack of nestedness of the 
discrete spaces. Convergence then follows from a weak reliability estimate of Qt, namely 

(3) Qe->0 => \\u - U t \\Hi(fl) 0, 

since qh < A^ — > as £ — > oo. We point out that our convergence proof makes use of the 
so called estimator reduction and does therefore not need the interior node property, which 
makes the result fairly independent of the local mesh-refinement strategy. 

1.3. Outline of current work. In Section 2, we formulate the continuous model problem 
and recall its unique solvability. In Section 3, the same is done for the discretized problem. 
Section 4 is a collection of the main results of this paper. Here, we introduce the error 
estimator g#, which is a generalization of the corresponding estimators from [9, 32]. We then 
state its weak reliability (Theorem 6) and our version of the adaptive algorithm (Algorithm 
7). Finally (Theorem 8), we state that the sequence of discrete solutions indeed converges 
towards the continuous solution u G H 1 ^). The subsequent Sections 5-7 are then devoted 
to the proofs of the aforementioned results and numerical illustrations. 

2. Model Problem 

2.1. Problem formulation. We consider an elliptic obstacle problem in M. 2 on a bounded 
domain Q with polygonal boundary T := dQ. An obstacle on Q is defined by the smooth 
function \ H 2 (Vl). Moreover, we consider inhomogeneous Dirichlet data g G H l l 2 (T) and 
thus additionally require \ < g almost everywhere on T. By 

(4) A := {v G H 1 ^) : v > x a.e. in f2, v\r = g a.e. on T}, 

we denote the set of admissible functions. For given / G L 2 (Q), we consider the energy 
functional 

(5) J(v) = l(v,v))-(f,v), 
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with the bilinear form 

(6) (u , v)) = / Vu-Vvdx for all u,v e H\Q) 
and with the L 2 -scalar product 

(7) (f,v)= [ fvdx. 

Jn 

By I ■ I, we denote the energy norm on Hq(Q) induced by ((■ , •)). The obstacle problem then 
reads as follows: Find u G A such that 

(8) J(u)=mmJ(v). 

VGA 

2.2. Unique solvability. For the sake of completeness and to collect the main arguments 
also needed below, we recall the proof that the obstacle problem (8) admits a unique solution. 
We stress that the following argument holds for any measurable obstacle \ with meaningful 
trace x|r- Our restriction to smooth obstacles x H 2 (Q) is needed for the equivalent 
reformulation in Section 2.3. 

Proposition 1. For given data (%, g, /) G H 2 (£l) x H l / 2 (Y) x L 2 {Vt), the obstacle prob- 
lem (8) admits a unique solution u E A which is equivalently characterized by the variational 
inequality 

(9) {u,v-u} > (f,v-u) 

for all v G A. □ 

The following two lemmata provide the essential ingredients to prove Proposition 1. We 
start with a well-known abstract result, cf. e.g. [6, Section 2.4-2.6]. 

Lemma 2. Let H be a Hilbert space with scalar product ((• , ■)) and IC C "H be a closed, 
convex, and non-empty subset. Then, for given L G %* , the variational problem 

(10) J{u) = min J(v) with J(v) = -({v , v) - (L , v) 

has a unique solution u G K,, where (■ , ■) denotes the dual pairing between % and %* . In 
addition, this solution is equivalently characterized by the variational inequality 

(11) (u, v-u))>(L,v-u) 

for all v G K. □ 

To apply Lemma 2, we observe that the obstacle problem (8) can be shifted into a setting 
with homogeneous Dirichlet data and % = Hq(Q). This involves a standard lifting operator 
C : H l l 2 (T) — > if 1 (f2), see e.g. [26, Theorem 3.37], with the properties 

(12) (£v)\ r = v and ||£u||.ffi(n) < Ci ||u||^i/2(r), 

where the constant C\ > depends only on Q. 

With elementary algebraic manipulations, see e.g. [24, Section II. 6], one obtains the fol- 
lowing well-known link between the obstacle problem (8) and the abstract minimization 
problem (10) from Lemma 2. 

Lemma 3. Let g G H 1 ^) be an arbitrary extension of g E H l l 2 {Y), e.g., ~g — Cg. For u G A 
and u — 'g — u G K, :— {v G Hq(Q) : v > x ~ d}> the following statements are equivalent: 
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(i) u solves the obstacle problem (8). 

(ii) u satisfies ((u , v — u} > (/, v — u) for all v G A. 

(iii) u satisfies ((« , v — u} > (/, v — u) — ((g , v — u} for allv G K,. □ 



2.3. Reduction to problem with zero obstacle. The following proposition allows 
to restrict to obstacle problems (8), where the obstacle x is even zero. This provides the 
formulation which can be treated by our adaptive method below. 

Proposition 4. For some smooth obstacle x £ H 2 (VL), the obstacle problem (8) with data 
{Xi9if) i s equivalent to the obstacle problem with data (0,g — x\r, f + If u G H 1 ^) 

solves the obstacle problem with data (0,g — x|r, / + u + x is the unique solution with 
respect to the data (x, g,f)- 

Proof. The solution u G A of the obstacle problem with data (x, g, f) is characterized by 

(u , v - uj > (/, v - u) 

for all v G A — {v G -ff 1 : v > x a - e - m ^"^Ir = g}- Substitution u := u — x an d 
v := v — x shows that this is equivalent to 

((u , 5; — S)) > (/, u — w) — ((x , u — u)) for all ? e i, 

where i := {? 6 : u > a.e. in Q,v\r = g — x\r}- Finally, integration by parts with 

w := v — u G (fi) proves 

— ((% , u — m}} = — / V% • Vw dx = Ax w dx. 

Jn Jn 

This concludes the proof. □ 



2.4. Model problem. According to the observation of Proposition 4, we may restrict 
to the case x = in the following. Having obtained a FE approximation Ue of u for the 
zero obstacle case, we may simply consider Ue, + x to obtain an approximation of the original 
problem with obstacle x £ H 2 {Vt). Since our analysis directly covers affine obstacles, we shall 
allow that x is globally affine on Q, i.e. we consider Problem (8) with respect to the data 

fcjJlGP'fflix^rjx^fl). 



3. Galerkin Discretization 



3.1. Problem formulation. For the numerical solution of (8) by an adaptive finite 
element method, we consider conforming and in the sense of Ciarlet regular triangulations 
Ti of Q and denote the standard Pl-FEM space of globally continuous and piecewise affine 
functions by l S 1 (7f). Note that a discrete function Ve G S l {%) cannot satisfy the inho- 
mogeneous Dirichlet conditions g G H l l 2 {Y) in general. We therefore have to approximate 
g ~ gi G 5 1 (7f|r), where the space 5 1 (7£|r) denotes the space of globally continuous and 
piecewise affine functions on the boundary V. For this discretization, we assume additional 
regularity g G -f^ 1 (r) and consider the approximation gi G S l (Te\r) which is derived by nodal 
interpolation of the boundary data. Note that nodal interpolation is well-defined since the 
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Sobolev inequality on the ID manifold T predicts the continuous inclusion H l (T) C C(T). 
Altogether, the set of discrete admissible functions At is given by 

(13) At := {Vi G S\%) : V e > in fi, V e = gt on T}, 
and the discrete minimization problem reads: Find Ut G Ae such that 

(14) J(U e ) = min J{V t ). 

Remark. Note that < g G H 1 ^) C C(r) implies that < (7(2;) = gf^(z) /or a// nodes 
z G T. Therefore, we conclude < ge on V for the nodal interpolant ge G 5 1 (7^|r)- D 



3.2. Notation. From now on, A/^ denotes the set of nodes of the regular triangulation 
%. The set of all interior edges E = T + n T~ for certain elements T + ,T~ G 7^ is denoted by 
£{\ The set of all edges of % is denoted by Et- In particular, Ee\v := := £e\£f contains 
all boundary edges and provides some partition of V. 

We recall that i5 1 (7^) = \Vn G C(f2) : Vi\t affine for all T G 7^} denotes the conforming 
Pi-finite element space. Moreover, S^(T e ) = S 1 {%) n Fq 1 ^) = {V e e S l {T e ) : V^|r = 0}. 

3.3. Unique solvability. In this section, we recall that the discrete obstacle problem (14) 
admits a unique solution which is again characterized by a variational inequality. 

Proposition 5. The discrete obstacle problem (14) admits a unique solution U# G Ai, which 
is equivalently characterized by the variational inequality 

(15) {U t , V t -Ui} > (f, Vt-Ut) 

for all V e G Ai. □ 

The proof of Proposition 5 is obtained as in the continuous case. It relies on the fact 
that discrete boundary data gi G iS 1 (7f|r) can be lifted to a discrete function <^ G S 1 ^) 
with <%|r = gt- The proof of the latter is a consequence of (and even equivalent to, see [17]) 
the existence of the Scott-Zhang quasi-interpolation operator P e : H 1 (Q) — > S 1 ^), see [35], 
which is a linear and continuous projection onto iS 1 (7f), i.e. 

(16) Pl = Pi and \\Piv\\ H i {a) < C 2 \\v\\ HHa) , for all v G H 1 ^), 
that preserves discrete boundary conditions, i.e. 

(17) (P^)|r = 'u|r for all v G /J 1 ^) with v\ T E S 1 (%\ r). 

The constant C2 > depends only on the shape regularity constant cr(7^) and on the diameter 
of Q. Moreover, Pi has a local first-order approximation property which is, however, not used 
throughout. 

4. Adaptive Mesh-Refining Algorithm and Main Results 

4.1. Newest vertex bisection. For the mesh-refinement, we use newest-vertex bisection. 
Assume that Mi C Ei is a set of edges which have to be refined. The refinement rules are 
shown in Figure 1, and the reader is also referred to [37, Chapter 4]. We stress a certain 
decay of the mesh-widths: 

• Marked edges E G A4e are split into two edges E', E" G £i + \ of half length. 
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Figure 1. For each triangle T G T, there is one fixed reference edge, indicated 
by the double line (left, top). Refinement of T is done by bisecting the reference 
edge, where its midpoint becomes a new node. The reference edges of the son 
triangles are opposite to this newest vertex (left, bottom). To avoid hanging 
nodes, one proceeds as follows: We assume that certain edges of T, but at least 
the reference edge, are marked for refinement (top). Using iterated newest 
vertex bisection, the element is then split into 2, 3, or 4 son triangles (bottom). 



• If at least one edge E of an element T G Ti is marked, T is refined into up to four son 
elements T G Te+i with area |T|/4 < \T'\ < |T|/2, cf. Figure 1. 
Moreover, given an initial mesh 7o, newest vertex bisection only leads to at most 4 • ^7o 
similarity classes of triangles. In particular, the generated meshes are uniformly shape regular 

diam(T) 2 

(18) sup a(Te) < oo with <J\Ti) = max — . 

leN Te7 « \T\ 

Furthermore, the number of different shapes of e.g. node patches that can occur, is finite. 
These observations will be necessary in the scaling arguments below. 

4.2. Weakly reliable error estimator. Let u G A denote the continuous solution of (8) 
and Ui G At be the discrete solution of (14) for some fixed triangulation %. To steer the 
adaptive mesh-refinement, we use some residual-based error estimator 

(19) g 2 := M^) 2 + °sc,(£) 2 ) + E ( a P x ^) 2 + osc,(£) 2 ) 

which has essentially been introduced in [9] for homogeneous Dirichlet data g = 0. First, 
r]i(E) 2 denotes the weighted L 2 -norms of the normal jump 

(20) Ve(E) 2 :=h E \\[d n U e ]\\ 2 LHE) hrEESf 

with He = diam(E') the length of E and [•] the jump over an interior edge E = T + f)T~ G E?. 
Second, for interior edges, oscg(E) 2 denotes the data oscillations of / 

(21) osc e (E) 2 :=\n i>E \\\f-f ne Jl 2M iovEeEf 

over the patch Q^e — T + U T~ associated with E, where the corresponding integral mean of 
/ is denoted by fn i E = (l/\£le tE \) J n f dx. Third, for boundary edges E G Ef and T G % 
the unique element with E C dT fl T, qi involves the weighted element residuals 

(22) osc,(£) 2 :=|T|||/||i 2(T) for E G Sf . 

Finally and in order to control the approximation of the nonhomogeneous Dirichlet data g 
by its nodal interpolant ge, we apply an idea from [3] and use 

(23) aV x e (E) 2 :=h E \\(g-g e y\\ 2 L2{E) for E G Ef, 
where (•)' denotes the arclength derivative. 
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To state a reliability result for the proposed error estimator g^, we need to introduce a 
continuous auxiliary problem. Given the discrete Dirichlet data ge G S 1 (Te\r), we define 

(24) A\ := {ve G H 1 ^) : v t > x a - e - i n ^> v e\r = Qt a.e. on T}. 
Applying Proposition 1 for the data (x,9i, f), we see that the auxiliary problem 

(25) J{u e ) = min J(v e ) 

admits a unique solution ue G A\. The following theorem now states weak reliability of ge in 
the sense that qi — > implies \\u — Ui\\ H i^ —> as I — > oo. The proof is given in Section 5 
below. 

Theorem 6. With Ue G A\ the solution of the auxiliary problem (25), the error estimator Q£ 
from (19) satisfies 

(26) i|K - us 2 < J(u e ) - J(u e ) < ciol 

with 

(27) qI := (Vi(E) 2 + osc £ (E) 2 ) + J2 osc e (E) 2 , 
as well as 

(28) ||u - Ui\\ H i {n ) < \\\u - u e \\\ + V2C 3 g e , 

where the constant C% > depends only on a (72) and on Q. Moreover, there holds 

(29) Qi^O =^ lltt-^H^^o. 

Remark. We stress that the reliability estimate (26) depends on the use of newest vertex 
bisection in the sense that only finitely many shapes of node patches and edge patches can 
occur. For red- green-blue refinement [37, Chapter 4], the equivalence of edge and node oscil- 
lations is open. The entire analysis, however, applies for a coarser error estimator, where 
edge oscillations are bounded in terms of the element residuals ||^£/||l 2 (t) ■ D 



4.3. Convergent adaptive mesh-refining algorithm. 

of the adaptive algorithm in the usual form: 



We can now state our version 



Solve 



Estimate 



Mark 



Refine 



Throughout, we assume that the Galerkin solution Ue G S 1 (Te) is computed exactly. For 
marking, we use the strategy proposed by Dorfler [18]. 

Algorithm 7. Fix an adaptivity parameter < 9 < 1, let % with I = be the initial 
triangulation, and fix a reference edge for each element T G %■ For each I = 0, 1, 2, . . . do: 

(i) Compute discrete solution Ue G At. 

(ii) Compute r]e(E) 2 , osce(E) 2 , and a,px e (E) 2 for all E G 

(iii) Determine (minimal) set Aie Q &t which satisfies 

(30) 6q 2 < Yl (vi(E) 2 + osc e (E) 2 )+ ^ (apx,(£) 2 + osc,(£) 2 ) . 

(iv) Mark all edges E G M.i and obtain new mesh 7~e+i by newest vertex bisection. 
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(v) Increase counter 1 1— > £ + 1 and iterate. □ 

The following theorem states our main convergence results. 

Theorem 8. Algorithm 7 guarantees the existence of constants < k, 7 < 1 and A > and 
a sequence ae > with lim = stzc/i that the combined error quantity 

£->oo 

(31) A, := JXE/*) - J(u e ) + 1Q ] + A apx| > 
satisfies contraction up to the zero sequence ae, i.e. 

(32) A m < nA e + a e for all £eN. 
In particular, this implies 

(33) = lim A e = lim g e = lim ||u - Ut\\jp-(n) 
as well as convergence of the energies 

(34) J(u) = lim J(u e ) = lim J(U e ). 

i— >oo I— >oo 

5. Proof of Theorem 6 (Weak Reliability of Error Estimator) 

5.1. Stability of continuous problem. For the finite element discretization, we have 
to replace the continuous Dirichlet data g G H l / 2 (T) by appropriate discrete functions ge. To 
make this procedure feasible, we have to prove that the solution of the obstacle problem (8) 
depends continuously on the given data. 

In the following, let % be a Hilbert space with scalar product ((• , •)) and /C, /Q C % be 
closed, convex, and non-empty subsets of %. We say, that /Q converges to K in the sense of 
Mosco if and only if 

(35) Vi> G /C 3v£ G /Q : V£ — > v strongly in "H as £ — > 00 
and 

(36) Vi> G "H Wvi G /Q : (f^ — i> weakly in % as £ — > 00 =>- v G /C) . 
We can now state an abstract stability result of Mosco: 

Lemma 9 ([27, Theorem A]). Assume the sets /C, /Q satisfy (35) -(36). Let L,Lp G H* be 
functional and assume Li — > L m "H* as £ — >■ 00. Finally, let u G /C and G /Q &e the 
unique solutions of (10) with respect to the data (JC, L) and (/Q, Lg), respectively. Then there 
holds strong convergence 

(37) U£ — >■ u in % as £ — >■ 00, 



ie. £/ie solution of the variational problem (10) depends continuously on the given data. □ 
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Remark. We remark that we have stated [27, Theorem A] only in a simplified form. In 
general, the preceding lemma of Mosco includes the approximation of the bilinear form ((• , ■)) 
as well, and it also holds for nonlinear variational inequalities, where the underlying operators 
Ai,A : H —> H* , e.g. Av := ((i> , •)) , are assumed to be Lipschitz continuous and strictly 
monotone with constants independent of t. □ 

We are now in the position to show that the solution u G A of the obstacle problem (8) 
continuously depends on the boundary data g G H 1 ^ 2 ^). 

Proposition 10. Recall that ue G A* t denotes the solution of the auxiliary problem (25). 
Provided convergence gi — > g in H 1 ^ 2 ^) of the Dirichlet data, there also holds convergence 
Ui — >■ u in H l (VL) of the (continuous) solutions as £ — > oo. 

Proof. As in the proof of Proposition 1, we define the sets 

K = {v G H%(n) : v> X -Cg} and K,* £ = {v e G H^Sl) : v e > X - Cgt}. 

We stress that both sets are convex, closed, and non-empty subsets of Hq(Q). As usual, let 
:= Hq(Q)* denote the dual space of Hq(Q). With respect to the abstract notation, 
we have L, Le G H~ 1 (Q) with 

(L , v) = (/, v) - (Cg , v)) and (L t , v) = (/, v) - (Cg e , v} . 

Note that this implies 

\\L - Li\\ H -na) < \\Cg - CgtWni^p) < \\g - 5^||ifi/2(r) ^°°> 0. 

Next, we prove that K\ — > JC in the sense of Mosco (35)-(36). 
To verify (35), let v G K, be arbitrary. Define 

v := v + Cg G A 

and let we G H 1 (Q) be the unique solution of the (linear) auxiliary problem 
Avj£ = in Q with Dirichlet conditions wg = ge — g on V. 
In the next step, we consider 

v e := v + w e G H 1 ^). 
Then, there holds Vi\r — ge as well as 

II — II II II ^11 II <?— >oo „ 

\\v - v e \\ H i in) = \\w £ \\ h i(q) < \\g e - g\\ H ^(T) > 0- 

Since the maximum of if 1 -functions belongs to H 1 , we obtain 

Vi := max{x,^} G H X (Q). 

By definition, there holds v e > x almost everywhere in Q as well as v#\r = gt since vi\v = 
9e — x\t almost everywhere on T. Altogether, this proves V£ G At- Moreover, since the 
function max{x, •} : H 1 ^) — > i/ 1 (f2) is well-defined and continuous (see e.g. [21, Chapter 
II, Corollary 2.1]), there holds 

v = max{x, v} = lim max{x, ve] = lim vg in iJ 1 (fi). 



Finally, we observe 

v e := v e 



Lgi e ~^°°) v — Cg = v in if 1 (fi) 
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and vi G /Q. 

To verify (36), let vi G K\ and assume that the weak limit vg — v v G -ffo(^) exists as 
£ — > oo. According to the Rellich compactness theorem, there is a subsequence {v£ k ) which 
converges strongly to v in L 2 (Q). According to the Weyl theorem, we may thus extract a 
further subsequence (vg k ) which converges to v pointwise almost everywhere in Q. Moreover, 

by continuity of the lifting operator, there holds C>gt k . — > Cg in H l (Vt) as j — > oo. Since 

this implies L 2 -convergence and again according to the Weyl theorem, we may choose a 
subsequence {C>gt h . ) which converges to Cg pointwise almost everywhere in Q. Altogether, 
the estimate 

ve k . > X ~ £gi k . a.e. in Q 
and monotonicity of the pointwise limit imply 

v > x ~ £>g a.e. in Q, 

whence v G K. 

Now, we may apply Lemma 9 to see that the unique solutions u G /C and ui G /Q of 
the minimization problems (10) with respect to the data (/C, L) and (/Q,Z^), respectively, 
guarantee 

\\u - Ui\\ H i(Q) > 0. 

With Lemma 3, there holds u = u + Cg as well as ii£ = U£ + Cgg. Altogether, we thus obtain 

I, || , |i~ — || || || l— >oo _ 

\\u - Ui\\ H i(p.) < \\u - ut\\ H i{^) + \\g - geWm/^r) > 

and conclude the proof. □ 



5.2. Some elementary preliminaries. The next result is an elementary lemma for 
variational inequalities, see e.g. [32, Proposition 2]. 

Lemma 11. Suppose that A is a convex subset of H l (VL) and f G L 2 (Q). Let U G Abe the 
solution of the variational inequality 

(38) (U,V- U)) >(f,V- U) for all V G A. 
Then, for all W G A, there holds 

(39) \\W ~W\f <J(W)-J(U), 

i.e. the difference in the H 1 -seminorm is controlled in terms of the energy. □ 

5.3. Proof of Theorem 6. We start with estimate (26), where the lower estimate clearly 
holds due to Lemma 11. We thus only need to show the upper estimate. To that end, we 
define a e G H~ l (Q) by 

{a e , v) := (f, v) - (u e , v)) for all v G #o(fi). 
Since ug\r = £/g|r, we ma Y argue as in the proof of [9, Theorem 1] to see 

UK - u 4\ 2 + ( a e 7 u e ~ U e ) < gj, 
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where the hidden constant depends only on the shape regularity constant cr(7}). A direct 
calculation finally shows 

J(U t ) - J(u £ ) = h$Jt , U e } - (f, Ui) - {^(u e , u £ } - (f, utj) 

+ ((ue , u e - U e } - (/, ue - Ue) + (a e , u e - Ue)) 
= -\\\ui - U e \\\ 2 + {<r t , u e - U £ ), 

whence J{U e ) - J(u t ) < gf. 

Next we prove (28). According to the Rellich compactness theorem and the triangle in- 
equality, there holds 

\\u - Ue\\m(o.) — t u ~ Uif + \\g - gi\\ H v*{T) 

< \\\u - Ui\\\ + \\\u£ - Ue\\\ + \\g - ge\\m/2(r)- 

From (26), we infer 

IN - US 2 ^ Qt < 
From [3, Lemma 3.3], we additionally infer that 

\\g - geWm/^r) < a P x ^( E ) 2 ^ e?» 

and the constant depends only on V and the local mesh ratio of Sf = 7}\r, whence on Q and 
<j(7~e). The combination of the last three estimates yields (28). Moreover, according to the 
last estimate, the convergence Qe — > implies gt — > g in H l l 2 (T) as £ — > oo. According to the 
stability result from Proposition 10, this yields convergence ue — > u in H l (Q), and we also 
conclude the proof of (29). □ 

6. Proof of Theorem 8 (Convergence of AFEM) 

The idea of our convergence proof is roughly sketched as follows: Dorfler marking (30) 
yields that Qe is contractive up to \\Ue+i — E^|||, see Section 6.1. Since there seems to be 
no estimate for this term to be available (recall nonconformity At % -Ae+i), we introduce a 
discrete auxiliary problem with solution Ue+i/ which allows to control 

|||^+l ~ Ue\\\ < lll^+l ~~ f^+l/lll + |||^M-l,<? ~~ ^111 

in terms of the energy J{Ui) — JiUe+i^e), see Section 6.3. In Section 6.4, we finally follow 
the concept of [15] to prove Theorem 8. 

6.1. Estimator reduction. The following contraction estimate is called estimator re- 
duction in [2, 15]. 

Proposition 12. Suppose that the set M.e Q &e satisfies (30) and that marked edges are 
refined as stated in Section J^.l. Then, there holds 

(40) Qj +1 <qej + C 4 \\\U i+1 -Ue\\\ 2 

with some contraction constant q G (0,1) which depends only on 6 £ (0,1). The constant 
C*4 > additionally depends only on the initial mesh 7q. 
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Since the proof of Proposition 12 follows along the lines of the proof of [32, Proposition 3], 
we only sketch it for brevity. 

Sketch of proof . First, the Young inequality proves for arbitrary 5 > 

Ql +1 < Y °sc m (£') 2 + Y a P x m(^') 2 
E'e£i+i E'e£f +1 

(41) 

+ (1 + 5) h E , \\[d n Uz}\\ 2 L2{EI) + (1 + 8- 1 ) h E , \\id n (U e+1 - U e )]\\ 2 LHE , y 

E'eep +1 E'e£f +1 

A scaling argument allows to estimate the last term by 

(42) Y h E/ \\[d n (U i+l -Ui)]\\l 2(EI) <C\\\U i+1 -U4 2 , 

and the constant C > depends only on o~(J~i). Next, one investigates the reduction of the 
other three terms if the mesh is refined locally: According to [32, Lemma 5], it holds that 

(43) Y h * WI^UMUe') < Y ^(Ef -\ Y ^(E) 2 . 



-t+i 



Next, [32, Lemma 6] resp. [31, Lemma 3.9.6] yields 

(44) Y osc m 

osci(E) . 

E'e£ i+1 E££ t Ee£ e nM t 

Moreover, it is part of the proof of [3, Theorem 4.2] that 

(45) Y apx m (£') 2 < Y ^{E) 2 -\ Y We( E ) 2 . 

E'e£f +1 Eeef Ee£fnM e 

Combining the estimates (41)-(45) and using the Dorfler marking (30), we easily obtain 

£? 2 +i <(1 + 6)(1- 6/4) g 2 + (1 + 5~ X )C fU t+1 - U t f. 

Finally, we may choose S > sufficiently small to guarantee q := (1 + 5)(1 — 6/4) < 1 to end 
up with (40). □ 

Remark. In the conforming case Ai C Ai + \, it is easily seen that the sequence Ui of discrete 
solutions tends to some limit which is the unique Galerkin solution with respect to the 
closure of {J^ Ai, see [31, Lemma 3.3.8]. Therefore, lim^ |||E^+i — L^||| = ; and elementary 
calculus predicts that (40) already implies lim^ qi = 0, cf. e.g. [2, Proposition 1.2]. Then, 
Theorem 6 would conclude thatlime \\u — UiWh 1 ^) = 0, i.e. u = U^. — Now that An % At + \, 
we did neither succeed to prove that the a priori limit Uoo exists nor the much weaker claim 
that \im.£ \\\Ui + i — L^||| = 0. □ 



6.2. Some a priori convergence results. In the adaptive algorithm, the mesh 7^+i is 
obtained by local refinement of elements in %. Consequently, the discrete spaces S l (%) are 
nested, i.e., 

(46) S l (Tt)QS 1 {Ti+i) foralUeN 



CONVERGENCE OF AFEM FOR OBSTACLE PROBLEM 13 

and the analogous inclusion also holds for the spaces 5 1 (7^|r) C ^C^+ilr) on the boundary. 
The following lemma is part of the proof of [1, Proposition 10]. 

Lemma 13. The nodal interpolants g% G <S 1 (7^|r) of some Dirichlet data g G i? 1 (r) converge 
to some a priori limit in if 1 (r), i.e. there holds 

/ A rr\ II II £—>OQ ~ 

(47) Hs-oo - gi\\H\T) > 

for a certain element g^ G H l (T). □ 

In the following, we will only use the convergence of ge to g^ in H l l 2 (T) as well as the 
uniform boundedness sup ||s^||ffi/2(r) < oo which is an immediate consequence. 

Proposition 14. The sequence ui G A\ of solutions of the continuous auxiliary problem (25) 
converges to some a priori limit Uoo G H 1 ^), i.e. 

(48) 1 1 it,,,, - ut\\w(si) — ^ 
for some function G H 1 ^). In particular, there holds 

(49) J( Uoo ) = lim J(u t ). 

i— >oo 

Proof. Since convergence in H l ' 2 (T) implies convergence in L 2 (T), the Weyl theorem applies 
and proves that a subsequence (gt k ) converges to g^ pointwise almost everywhere on T. 
Therefore, the limit function g^ G iJ 1 ^ 2 (r) satisfies g^ > x almost everywhere on T. In 
particular, the obstacle problem (8) with g = g^ has also a unique solution G H l (T). 
Applying the stability result from Proposition 10, we obtain convergence of ue to in if 1 (f2) 
as £ — > oo. □ 



6.3. Discrete energy estimates. The main difficulty which we suffered in our analysis 
is that it is not clear that |||o^ +1 — Ue\\\ resp. J(Ui + \) — J{Ui) tend to zero as t — > oo. We 
circumvent this question by introducing a discrete auxiliary problem: Find U^+ij G Ae + i/ 
such that 

(50) J(U t+ltt ) = ¥/ min J(V £+1 ). 
where the admissible set reads 

(51) Ai+ij := {V i+1 G S 1 (%+i) : V e+ i > x in Ve+i = geonT}. 

Applying Proposition 5 to the data (Ae+i^ge) instead of (Ae,ge), we see that the auxiliary 
problem (50) admits a unique solution Ui+ij G Ae+x s £. 

Our first lemma is a key ingredient of the upcoming proofs. It states that one may change 
the boundary data of a discrete function and control the influence within Q. 

Lemma 15. Let Wg + \ G S 1 (% +1 ) with W i+1 \r = gi+x- Define Wf +1 G S 1 (% + i) by 



(52) WUz) 



W e+1 (z) forzeAf e+1 \T, 
ge(z) for z G Af i+1 n T. 

Then, with the local mesh-width hi G L°°(r) defined by H(\e = h,E, there holds 

(53) \\W e+1 - We +l \\ H i {n) < C 5 Wh 1 / 2 (g e+l - gi)'\\ L *{T) < C 6 \\gt+i ~ 9i\\hV*(T)> 
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where C§,Cq > depend only on % and on Q. 

Proof. By definition, we have Wf +1 G 5 1 (7£+i) with Wf, Jr = 9e- Moreover, there holds 
W/+iIt = W^+i|t for all T G %+i \ T[ +l . Let T = convex, z 2 , z 3 } G T[ +l and $ T : T rcf -)> T 
be the affine diffeomorphism with reference element T rc f = conv{(0, 0), (0, 1), (1, 0)}. 

The transformation formula and norm equivalence on P 1 (T re f) prove 

\\W l+1 - W[ +1 \\l*(t) = |T| 1/2 \\{W i+1 - Wj +1 ) o $t||l»(2W) 

3 
5=1 

<|fi| 1/2 ^|(W m -M// +1 )fe)l 

as well as 

||V(Wi+i - V^ +1 )|| i2(T) < a(%) \\V((W e+1 - W l i+1 ) o $ T ) || i2(Tref) 



<a(7^)^|(W m -W/ +1 )0 

Now, we consider the nodes z G {zi, 22, ^3} of the triangle T: For z G A/« ft T holds <%(z) = 
5^+1 ( 2 ) by definition of the nodal interpolants. Therefore, 

(Wi+i - wtf +1 )(*) = if ^ g (jvi+i n n) u n r). 

Thus, it only remains to consider nodes z G (A^+i\A/^) D T. In this case, z is the midpoint 
of a refined edge of the unique father T G 7^ with T C T, and there even holds £ z cT. 
Let z G A/i be an arbitrary endpoint of i? 2 . Then, 

(W e+1 - Wt +1 )(z) = (g i+1 - g e )(z) = 0, 

and the fundamental theorem of calculus, now applied for the arclength derivative, yields 

\{W l+1 -W* +1 ){z)\< [jge +1 -ge)'dT <\z-^ 1/2 \\i9i-g^i)'\\^s a y 

Jz 

Combining our observations, we have now shown 

1 /2 

||^ +1 -^/ +1 b 1(T) < hf\\(gt-gt + i)'\\ L > {2 )-( E Mfo-ft+O'll^ ' 



-ES£i B6£i 



ECdf ECdT 



Finally, note that each edge E G Sf belongs only to one element T G %■ Thus, E can at 
most be selected by all (but at most four) sons T G Te+i of T, cf. Figure 1. This observation 
yields 

\\w e+ i-w* +1 \\ 2 HH n)= E II^ + i-^/ + iIIh 1( t) < E^ 2 !!^-^ 1 )'!!^)- 

With the local mesh-width function he G L°°(T), this estimate reads 

\\W l+ i - 11^-1(0) < \\h\ /2 (ge - <fc+i)'IU a (r) < ||& - 9e+i Wh 1 / 2 ^)^ 
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where we have used a local inverse estimate from [14, Proposition 3.1] in the last step. We 
stress that the constant depends only on the local mesh-ratio of Te+i\r, whence on a(Te+i). 
This concludes the proof. □ 

With the notation of Lemma 15, we can now formulate an additional convergence result. 
Proposition 16. The sequence of discrete solutions Ue + \/ G At+i/ satisfies 

(54) \J(U W ) - J(U t+1 )\ + \J(Ui +lte ) - J(U e e+1 )\ < C 7 \\g e+1 - g e \\ H i/ m ^ 
with some constant C-j > which only depends on Tq. 

To prove Proposition 16, we start with the observation that the energy functional J7"(-) is 
coercive. 

Lemma 17. For each sequence uy G if x (f2) with changing boundary data wi\r = gt, bounded 
energy implies uniform boundedness in H 1 ^), i.e. 

(55) sup Jiwz) < oo =>- sup 1 1 iu^I | #i < oo. 

Proof. Note that triangle inequalities and the Friedrichs inequality yield 

||^IU 2 (n) < \\we - £>9e\\L 2 {n) + II^IU^o) 

< \\V(wi - Cg e )\\ L 2 (n) + \\£gi\\ L 2 {n) 

< ||V^||i2(Q) + \\Cge\\H l {n) 

< \\Vw e \\ij2(n) + M 

with M = sup e WgiWH 1 / 2 ^) < 00 > where we have finally used continuity of £ and a priori 
convergence of ge, cf. Lemma 13. Consequently, we obtain 

IMIifi(fi) < ||Vu;/||i3(n) + M 

and therefore 

JM = -||V^||| 2(n) - (f,w t ) > (||^||hi(o) - M? - ||/|U2 (n) ||^|| H i (n) . 
The implication (55) is an immediate consequence. □ 

Lemma 18. The sequence of discrete solutions Ut G At from Algorithm 7 satisfies 

(56) sup \,J{Ut)\ < oo as well as sup < oo. 

Proof. According to Lemma 17 with wi = Ue, it is sufficient to prove boundedness of the 
energy. To that end, we define 

V e := max e {P e £g<> - X, 0} + x G Ai, 

where maxj denotes the nodewise max-function, i.e. maxi{v,w} G 5 1 (7^) is defined by 

maxg{v, w}(z) := max{t> (z), w(z)} for all z G Afe- 

We consider the function We := ma.xe{PeCge — X;0}. Note that, by definition, < 
\(Pe£ge — x)( z )\ f° r a U z e -A/i- According to a standard scaling argument and PeCge — x £ 
S 1 ^), we infer, for all T G %, 

\\W e \\ L 2 {T) < \\P t Cg t - xlU 2 (T) as well as || VW e \\ L 2 {T) < \\V(P e Cg e - x)IU 2 (t) 



16 M. FEISCHL, M. PAGE, AND D. PRAETORIUS 

due to the nodewise estimate. From this, we obtain a constant C > with 

< \\Pg£ge - xll-ff 1 ^) + llxll// 1 ^) ~ IMI#V2(r) + Hxll h 1 ^) — C < 00 
since sup £ ||<%||ffi/2(r) < oo and the operator norm of Pi depends only on cr(7~e). Therefore, 

J{Ut) < J(V e ) = l - ||VF,||i 2(n) - (/, Vt) <C(C+ ||/|U 2(n) ) < oo 
yields sup J{Ui) < oo and concludes the proof. □ 



Proof of Proposition 16. According to Lemma 18, we have M := sup £ < oo. 

With the aid of Lemma 15, we obtain 

(57) 

J(U e+1/ ) < J(U e w ) = i \\\u e e+1 \\\ 2 ~ (/, Ut-i) 

< l - (|||f/ m ||| + \\\U £ e+1 - U t+l \\) 2 - (f,U e+1 ) + ||/|| L2( n)||^ +1 - U i+1 \\» m 

= J(U i+1 ) + -|||^/ + i-f^+i||| 2 + |||^+illl III Ui +1 -U i+1 III + ||/|U2(n)||t/'l + i-C^+i||i2(n) 

< J(U l+ i) + - C\ \\g t - ^+i||^-i/2 (r) + C 6 (M + ||/||x,2(n)) \\gi+x ~ giWn 1 / 2 ^)^ 
whence 

J(Ue +lt i) — J{Ui + i) < \\gi — g £+1 \\ 2 Hl/2 ^ + \\g i+1 — gi\\H^(r)- 

Estimate (57) now reveals 8\ip e J'(Ui+i t i) < oo, whence sup £ < °o according to 

Lemma 17 with W£ = Ue + \^. Arguing as in (57) with the ansatz J{Ui+\) < JiU^lf) shows 

J{Ui+\) — J(U e+ i : i) < \\gi — ge+i\\ 2 H i/2^ + \\ge+i — geWm/^fv)- 
The combination of the last two inequalities yields 

\J(U i+1 ) - J(U e+1 ^)\ < \\gt - ^+ill^i/2(r) + \\ge+i ~ 9e\\m/*(r) °- 

Finally, we recall J{Ui+\,i) < J{U^ +1 ). Then, the above calculation also yields 

\J(U e i+1 ) - J{U W )\ = J{Ul x ) - J{U W ) 

= [J(Uf +l ) - J(U e+1 )] + [j(U e+1 ) - J{U W )} 

<~ \\ II 2 i II II l-^oo „ 

^ \\9i — fl , «+i|liji/2( r ) + \\9i+i — ^||i?i/2(r) > u - 

This concludes the proof. □ 



We are now ready to prove that the sequence Ut G At of discrete solutions generated by 
Algorithm 7 indeed converges towards the exact solution u G A. 
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6.4. Proof of Theorem 8. With the help of Lemma 15 and Lemma 11 applied twice 
for Ug G Ai + i t £, and U$ +1 G Ai+i,e, we obtain 

2 lll^+i ~ UilW 2 < \\\Ui + i — Ue+i^Wl 2 + |||L^ + i 5 £ — Ui\\\ 2 

< 2 \\\U t+ i - Uf +1 f + 2 \\\U{ +1 - U w \f + 2 [J(^) - J(f/, +M )] 

< 2C 2 \\h 1 / 2 (g i+1 - g t )'\\h { v) + 4 [J(Uf +1 ) ~ .70W] 

+ 2[J(Ut)-J(U w )]. 

Recall that Ag + i = ^J(Ue+i) — ^(i^+i)] + 7^ +1 + Aapxf +1 . We now use the last estimate 
to see 

[j(U e+1 ) - J(u e+1 )] = [j(U t ) - J(u e )] + [J(u e ) - J(u e+1 )] 

+ [j(U e+1 ) - J(U t+ltt )] + [J(U £+1/ ) - J(U t )] 
< [j(U t ) - J(u e )] + [J(u e ) - J(u £+1 )] 

+ [j(U t+1 ) - J(U W )] + C 2 \\h\' 2 {g l+l - gdVmv) 

- i fU i+1 - U e \\\ 2 + 2 - ^(C^+m)]- 

Note that according to Proposition 14 and Proposition 16 it holds that 

ai := \j(u t ) - J{u i+1 )\ + \j(U t+1 ) - J(U l+1/ )\ + 2 - J(£W)| 

tends to zero as I — > oo. Next, we recall that on the ID manifold T, the derivative g' e of 
the nodal interpolant is the elementwise best approximation of the derivative g' by piecewise 
constants, i.e. 

(58) \\h) /2 (g - ^+i) / ||l 2( r) + \\h\ ,2 {g t+x - 9i)'f L \T) = \\ h \ /2 (9 ~ 9i)'\\h{v) 

according to the elementwise Pythagoras theorem. So far and with A = C 2 , we thus have 
derived 

A m < [j(Ui) - J(u e )] + 7^ +1 + Aapx^ +1 + A ||/iJ /2 (^ + i - gt)'\\h(r) 

- -\\\Ui +1 - U t \\\ 2 + a t , 

< [J(U e ) - J{ui)] + 7^+i + Aapx 2 , - - fU l+1 - U t \f + a t , 

where we have used the Pythagoras theorem in the second step. 

Next, the estimator reduction of Proposition 12 applies and provides constants < q < 1 
and C4 > with 

A m < [j(U t ) - J(ut)] + -fqg 2 e + Aapx 2 + ( 7 l7 4 - -) \\\U i+1 - [7 £ ||| 2 + a e 



< [J{U t ) - J(u t )] + jqg 2 + Aapx 



provided the constant < 7 < 1 is chosen sufficiently small. 
Next, from Theorem 6, we infer C^ 2 \ j{Ug) — J{ut)\ < g 2 . 
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We plug this estimate into the last one and use the fact that apx 2 < gj to see 

A m < (1 - 7 £C 3 - 2 ) [j{U t ) - J{u e )] + 7 (g + 2e) ^ + (1 - e^Aapx 2 . + a, 

with k := max{l — ^sC^ 2 , q + 2 e, 1 — 75}. For < 2 e < 1 — q, we obtain < n < 1 and 
conclude the proof of the contraction property (32). 

The application of Lemma 11 for At C A\ yields 3{ui>) < J{U() and hence A^ > 0. 
Together with ote > and a.g —> 0, elementary calculus yields A^ — > as I — > 00, cf. [2, 
Proposition 1.2]. From < Qi < A^, we obtain lim^ ^ = 0, and the weak reliability of Qi 
stated in Theorem 6 concludes the proof of (33). □ 




Figure 2. Galerkin solution Uq on adaptively generated mesh Tq with N = 
4.159 elements for 9 = 0.8. 



7. Numerical Experiments 

We consider numerical examples, one of which has also been treated in [9]. The mesh in 
each step is adaptively generated by Algorithm 7. For the solution at each level, we use 
the primal-dual active set strategy from [22]. The numerical results for example 1 are quite 
similar to those in [9] . We stress, however, that our approach includes the adaptive resolution 
of the Dirichlet data and, contrary to [9] , the upcoming examples are thus covered by theory. 



7.1. Example 1. We consider the obstacle problem with constant obstacle x — 011 
the square Q := (—1.5, 1.5) 2 and a constant force / = —2. The Dirichlet boundary data 
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Figure 3. Numerical results for uniform and adaptive mesh- refinement with 
9 = 0.8, where yfel, Qi, and apx^ are plotted over the number of elements 



g G H l (T) are given by the trace of the exact solution 

Y - ln(r) — \, r > 1 



u :- 



2 ""V I 2 

0, else 



where r — \x\ and | • | denotes the Euclidean norm on M 2 . The solution is visualized in Figure 2. 
We compare uniform and adaptive mesh-refinement, where the adaptivity parameter 9 varies 
between 0.4 and 0.8. The convergence history for uniform and adaptive refinement with 
9 = 0.8 is plotted in Figure 3, where the error is given in the energy functional, i.e. 

(59) et=\j(U e )-J(u)\ 

and the Dirichlet data oscillations apx £ are defined by (23). Note that due to Theorem 8, 
the adaptive algorithm drives the error and thus also the energy ee to zero, whence it makes 
sense to plot these physically relevant terms. All quantities are plotted over the number of 
elements N = of the given triangulation. Due to high regularity of the exact solution, 
there are no significant benefits of adaptive refinement. We observe, however, that error 
and error estimator, as well as Dirichlet oscillations show optimal convergence behaviour 
0(N^ 2 ) and 0(N-^ 4 ) respectively. Also, the curves of Qa and y/JJ are parallel, which 
experimentally confirms classical reliability and efficiency of the underlying estimator qi in 
the sense of \\u — — Qi- 

Figure 4 compares different values of 9 and we can see that each choice of 9 G {0.4, 0.6, 0.8} 
leads to optimal convergence, i.e. the curves basically coincide. 
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Figure 5. Adaptively generated meshes 7s (left) and 7n (right) with A 
and A = 34.024 elements, respectively, for 9 = 0.6. 



678 



Finally, Figure 5 shows the adaptively generated meshes after 5 and 11 iterations. As 
expected, refinement basically takes place in the inactive zone, i.e. elements where the discrete 
solution Ug does not touch the obstacle. 
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Figure 6. Galerkin solution £7 10 on adaptively generated mesh 7io with N = 
8.870 elements for 9 = 0.6. 



10 1 




10 3 10 4 10 5 10 6 10 7 

Figure 7. Numerical results for uniform and adaptive mesh- refinement with 
6 = 0.4, where y/ej, Qi, and apx^ are plotted over the number of elements 
N = #%. 

7.2. Example 2. In this example, we consider the obstacle problem for a non- 
H 2 (Q) obstacle 

fi[sin(5(>+ (I-tt/IO))) + 1] x < -1, 
X I else. 
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10 3 10* 10 5 10 6 10 7 



Figure 8. Numerical results for yfel for uniform and adaptive mesh- 
refinement with 9 G {0.4,0.6,0.8}, plotted over the number of elements N = 




Figure 9. Adaptively generated meshes Tn (left) and 71s (right) with N = 
4.104 and N = 33.310 elements, respectively, for 9 = 0.6. 



on the L-shaped domain Q := (—2, 2) 2 \[— 2, 0). The right hand side is given in polar coordi- 
nates by 



f(r,(p) := -r 2/3 sin(2^/3)(d/dr(7i)(r)/r + d 2 /dr 2 j 1 (r)) 
A r -V3 d / dr ( 7l )( r ) s in(2^/3) - 72 (r), 
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where d/dr denotes the radial derivative. Moreover, r := 2(r — 1/4) and 

1, f<0, 
7i( r ) : = { -6r 5 + 15f 4 - 10f 3 + 1, 0<f<l, 
0, f>l, 



72 (0 := 



0, r < 5/4, 

1, else. 



The Dirichlet data g E H 1/2 (T) are given by the trace of the obstacle x- Since the exact 
solution for this problem is unknown, the Galerkin solution on a uniform mesh with approx- 
imately iV = = 1.500.000 elements has been used as reference solution. The non-affine 
obstacle was treated by means of Proposition 4. Again, we compare uniform and adaptive 
mesh-refinement for different adaptivity parameters 9 between 0.4 and 0.8. Figure 7 shows 
the convergence history for 9 = 0.4 plotted over the number of elements N = j^l~t. As before, 
adaptive refinement leads to the optimal convergence rates 0(N~ 1 ^ 2 ) and (9(iV~ 3//4 ) for ^fel 
and apx f , respectively. Due to the corner singularity of the exact solution at 0, we observe 
that uniform mesh-refinement leads to a suboptimal convergence behaviour. 

Figure 8 compares the error for uniform and adaptive mesh refinement, where the adapti- 
vity parameter 9 varies between 0.4 and 0.8. Again, we observe that each adaptive strategy 
leads to optimal convergence rates, whereas the convergence rate for uniform refinement is 
suboptimal. 

In Figure 9, the adaptively generated meshes after 13 and 18 iterations are visualized. As 
before, refinement is basically restricted to the inactive zone. 
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